Mon. Not. R. Astron. Soc. 000, 1-29 (2010) Printed 16th April 2010 (MN I^TeX style file v2.2) 



Solving the angular momentum problem in the cold 
feedback mechanism of cooling flows 

Fabio Pizzolato * 

INAF-Osservatorio Astronomico di Brera, Via Brera n.28, 20121 Milano (Italy) 

Noam Soker f 

Physics Department, Technion, Haifa 32000 (Israel) 
16th April 2010 



ABSTRACT 

We show that cold clumps in the intra-cluster medium (ICM) efficiently lose 
their angular momentum as they fall in, such that they can rapidly feed the 
central AGN and maintain a heating feedback process. Such cold clumps are 
predicted by the cold feedback model, a model for maintaining the ICM in 
cooling flows hot by a feedback process. The clumps very effectively lose their 
angular momentum in two channels: the drag force exerted by the ICM and 
the random collisions between clumps when they are close to the central black 
hole. We conclude that the angular momentum cannot prevent the accretion 
of the cold clumps, and the cold feedback mechanism is a viable model for a 
feedback mechanism in cooling flows. Cold feedback does not suffer from the 
severe problems of models that are based on the Bondi accretion. 

Key words: galaxies: clusters: intracluster medium galaxies: clusters: cooling 
flows galaxies: clusters: general galaxies: clusters: individual: Virgo 



1 INTRODUCTION 

For more than a decade now it is clear that the intra-cluster medium (ICM) in cooling 
flow (CF) clusters of galaxies and CF galaxies must be heated, and the heating process 
should be stabilized by a feedback mechanism (see reviews by Peterson & Fabian 2006 and 
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McNamara & Nulsen 2007). However, in many cases the heating cannot completely offset 
cooling (e.g., Wise et al. 2004; McNamara et al. 2004; Clarke et al. 2004; Hicks & Mushotzky 
2005; Bregman et al. 2006; Salome et al. 2008; Wilman et al. 2009) and some gas cools to 
low temperatures and flows inward (e.g., Peterson & Fabian 2006). The mass inflow rate is 
much below the one that would occur without heating, and the flow is termed a moderate 
cooling flow (Soker et al. 200f ; Soker & David 2003; Soker 2004). 

Two basic modes have been proposed to feed the super massive black hole (SMBH) at 
the heart of the active galactic nucleus (AGN) in CF clusters (e.g., Combes 2008). They 
can be termed the hot and the cold feedback mechanisms. In the hot-feedback mode the hot 
gas from the vicinity of the SMBH is accreted, such as in the Bondi accretion process (e.g., 
Omma & Binney 2004; Allen et al. 2006; Russell et al. 2010). More sophisticated schemes 
that use direct calculations of the inflow of the hot gas, but are basically similar to the Bondi 
accretion process, also exist (e.g., Ciotti et al. 2009, 2010). However, the Bondi accretion 
process suffers from two severe problems (Soker 2006; Soker et al. 2009), to the point that 
it fails. First, there is no time for the feedback to work. The accretion of gas from the very 
inner region does not have time to respond to gas cooling in the outer regions (r > 10 kpc). 
Second, although in galaxies the accretion rate can in principle account for the AGN power 
(Allen et al. 2006; Ciotti et al. 2009, 2010), in CF clusters the Bondi accretion rate is much 
too low (Rafferty et al. 2006). In a recent paper it was argued that the Bondi accretion 
process fails to maintain feedback in the process of galaxy formation as well (Soker 2009). 

In the cold feedback model (Pizzolato & Soker 2005; Soker 2006; Pizzolato 2007) the 
mass accreted by the central black hole originates in non-linear over-dense blobs of gas 
residing in an extended region of r ~ 5 — 30 kpc; these blobs are originally hot, but then 
cool faster than their environment and sink toward the centre (see also Revaz et al. 2008). 
The mass accretion rate by the central black hole is determined by the cooling time of 
the ICM, the entropy profile, and the presence of inhomogeneities (Soker 2006). In the cold 
feedback model if gas cools at large distances, so do the overdense blobs. In a relatively short 
time they feed the SMBH, and the ICM is heated on a time scale shorter than its cooling 
time scale. This overcomes the main problem of models based on the Bondi accretion process 
(Soker 2009). Wilman et al. (2009) suggest that the behaviour and properties of the cold 
clumps they observe in the cluster A 1664 support the cold feedback mechanism. In general, 
the presence of large quantities of cold (T < 10 4 K) gas in CF clusters (e.g., Edge et al. 
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2002; Edwards et al. 2007) suggests that the non-linear perturbations that are required to 
form the dense clumps do exist in CF clusters. 

The main arguments against the cold feedback mechanism was that the cold clumps 
supposed to feed the SMBH have too much angular momentum, and they cannot fall directly 
to the SMBH: for this reason, their accretion rate will be too slow and the response time 
too long (e.g., Russell et al. 2010). 

In this paper we show that dense clumps lose their angular momentum rapidly enough, 
such that there is no angular momentum problem in the cold feedback mechanism. The 
equations are presented in Sections 2 and 3, and the cluster properties used for the calcu- 
lations are given in Section 4. In Section 5 we present the solution for falling clumps, both 
without (Sec. 5.1) and with (Sec 5.2) angular momentum. In Section 6 the loss of angular 
momentum via collision is estimated. Our main conclusions are in Section 7. 



2 EQUATIONS OF MOTION 

Following Pizzolato & Soker (2005) (also Paper I), we write the set of equations governing 
the motion of a cold clump in the ICM. In Paper I we adopted the following equations for 
the motion of a clump: 

£ = v (la) 
p'V^ = -^fSp\v\v + g(p'-p)V, (lb) 

where the primed quantities refer to the clump. The first equation relates the position x of 
the clump with its velocity v; the second is Newton's second law for the clump's motion: V 
is the volume of the clump, and p' V its mass. The first term on the right-hand side is the 
drag force, where Cd ~ 0.75 (Kaiser 2003) is the drag coefficient and S the clump's cross 
section. The second term is the gravitational force, corrected for the buoyancy. 

These equations must be completed with the evolution of the density of the clump, which 
we give in a slightly different form than in Paper I. Let s(x, t) and s'(x, t) be the specific 
entropies of the ICM and of a clump at the time t and position x. We assume that the clump 
is always in pressure equilibrium with the ICM: P = P' . For this reason the specific entropy 
s' of an overdense clump is lower than the ICM specific entropy. Let 

As(x,t) = s'(x,t) -s(x,t) (2) 

be the entropy difference between the clump and the ambient ICM at the time t. We now 



4 Pizzolato & Soker 

consider the evolution of As at a later instant t + dt. On account of the radiative losses both 
the ICM and the clump have cooled. The entropy profile of the ambient gas at t + dt is 

s(x, t + dt)- s(x, t) = - Ue nH ^p T ' dt, (3) 

where A(T) is the cooling function, n e , nn and n being the electron, proton and total density 
of the ICM 1 . At the same instant t + dt the entropy content of the clump has changed as 
well due to the radiative losses: 

s'(x + dx, t + dt) - s'(x, t) = - W ° w jA(r') ^ (4) 

The position is now evaluated at x + dx because the clump is denser than the environment, 
so it has sunk a little, moving from x to x + dx. The entropy contrast with the ambient gas 
at the new position of the clump is 

As(x + dx, t + dt) = s'(x + dx,t + dt) - s(x + dx,t + dt). (5) 

Subtracting equations (5) and (2) taking into account equations (3) and (4) we find the 
variation of the entropy contrast: 

dAs = As(x + dx,t + dt) - As(x, t) = _ n e n ^( r dt + n en H A(T) ^ _ y ^ ^ 

n'T' nT 

The entropy contrast As may be expressed in term of the density contrast 

5 = (n' — n)/n (7) 



As = f£; B ln — - - | k B In — - = -| k B ln(l + 6); (8) 



under the isobaric assumption P = P' between the clump and the ICM: 

n /5/3 J 2 B "I I 5/3 J 2 

notice that if 5 = the entropy contrast vanishes. Equation (6) may then be written 

5 1 d< *_3 v vlnK , <n' H A(T') -n e n H K{T) ^ 



2/3 
e j 



2 l + 5dt 2 P 
where we have introduced the pressure P' = P = nk B T and the entropy index K = T/n 
connected to the specific entropy s through the relation 

s = lk B hiK. (10) 

The temperature of the clump is related to that of the ambient gas through the isobaric 
condition 

r = T /{i + §). (ii) 



1 In the present analysis we neglect the "background" flow due to the cooling of the ICM. The effect of this term was addressed 
by Loewenstein (1989). Its inclusion would require some modelling of the background flow. The impact on the results is weak, 
and does not justify the amount of required work. 
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Figure 1. Plot of the function lj(8, T) defined by Equation (13). 5 is the density contrast (eq. 7), while u is a function that 
enters the equation describing the relative evolution of the clump and the ICM; the larger ui is, the faster is the evolution of 
the clumps due to radiative cooling. Terms due to buoyancy and to radiative cooling appear in ui. The cooling function used to 
draw this plot has metal abundance Z = 0.3 (Anders & Grevesse 1989 solar units). The curves refer to different temperatures: 
T = 0.5 keV (solid line). T = 1.0 keV (dashed line) T = 2.0 keV (dotted line) and T = 4.0 keV (dash-dotted line). 



is plotted in Fig. 1 for several values of T. Note that for small values of 5 the profiles are 
very similar to each other, since 5 <C 1 implies co ~ (2 — din A/dlnT) 5 for any T. 

Since relaxed clusters are convectively stable, the entropy gradient is positive, and since 
v is directed inwards, the first term on the right-hand side of Equation (12) is negative, and 
tends to stabilise the clump, i.e. to keep its density contrast small. This term is responsible 
for the stability of clumps with small density contrasts. The second term describes the 
clump's instability to its own radiative cooling with respect to the ICM. It has the sign of 
co, which is usually positive for 5 > 0, thus giving a positive contribution to 5. As discussed 
in more detail in Appendix A, for small 5 this term is responsible for an overstability of the 
clump, as already found and discussed by Balbus & Soker (1989). 

As pointed out by Burkert & Lin (2000), when the cooling time of the clump falls below 
the sound crossing time, the pressure equilibrium does not hold any longer, and the clump 



With the aid of this constraint we may rewrite Equation (9) 




uo(6,T), 



(12) 



where /i# = n-n/n — 0.44. The function 




(13) 
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evolves isochorically. For consistency, we shall have to check a posteriori that the clump does 
not enter this regime during its evolution. 

Equation (12) is the same equation for the evolution of the density contrast 5 adop- 
ted in Paper I if the term with the ICM cooling is dropped. The new equation expresses 
more transparently than in Paper I the role of the ambient entropy gradient in the clumps' 
accretion (Soker 2006). The dynamical status of the overdensity (the right-hand side of 
Equation 12) depends on the properties of the clump only through its overdensity: the re- 
maining controlling parameters only depend on the ICM temperature, density and entropy 
distributions. At this stage, we are neglecting the ICM heating due to the AGN activity, and 
we only consider its radiative cooling. The ICM cooling term makes the clumps slightly more 
stable than considered in Paper I, where even a clump of zero overdensity gave a positive 
contribution to 5. 



3 THE ANGULAR MOMENTUM 

As in Paper I, we follow the evolution of a clump, but the purpose of this study is to address 
more fully the issue of the clump's angular momentum. If the clump is not able to lose 
effectively its angular momentum, it cannot be accreted by the AGN rapidly enough, and 
the feedback cannot work. 

We estimate the time scale taken by the drag force to dissipate the angular momentum 
of the clump. Taking the vector product of x by Equation (lb) we get 

- = -2 CdWI 1 (14) 
dt 8 a (1 + 5) { ' 

where a is the radius of the clump (assumed spherical) and 1 = x x v is the specific angular 
momentum. Note that since 1x1 = the orbit of a clump is planar. The characteristic time 

of angular momentum loss is 

8a(l + <5) 7 / a \ fl + 6\ ( o \~ l 

Tl = 1/1 = — A 1 ~ 1.2 x 10 7 yr —— I ) , 15 

1 3 C D |v| y ViOOpcy V 10 J VSOOkms- 1 / ' V ; 

where we have equalled |v| to a typical velocity dispersion a. This time is to be compared 
to the time of a Keplerian revolution to give 

Tl /r K ~(l + 8)a/r (16) 

For 5 ~ 10, a ~ 100 pc and r ~ 10 kpc this ratio is ti/tk — 10 _1 1. Therefore a clump of 
moderate overdensity is expected to quickly lose its angular momentum and fall straight on 
the AGN. An exception to this occurs when the cooling time of the clump is much less than 
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T\. In this case, the clump becomes very dense (8 ^> 1) before losing its angular momentum. 
For 5 ^> 1 the drag force is small, and the clump may orbit around the centre losing but 
slowly its angular momentum. This approximate derivation must be complemented by a 
more accurate numerical calculation, presented in Section 5. 

Since the orbit of a clump is planar, the differential equations for its (isobaric) evolution 
may be written in two-dimensional polar coordinates (r, <fi) 
dr 

— = u (17a) 
dt y ' 

£ = < 17b > 

I 2 s ^ l v l $ ,«„ \ 

"77 = ~7 ~ I —JZ ?T U — g- r (17c) 

dt r 3 8 a (1 + 5) 1 + 5 v ' 

s = -t c » sinij ' (17d) 



where u is the radial velocity, I the modulus of the angular momentum, g = |g| is the 
modulus of the gravitational acceleration, and 

| v | = (u 2 + l 2 /r 2 ) 1/2 (18) 

is the modulus of the velocity. 

The last equation we need governs the evolution of the clumps' size a. There are differ- 
ent possible choices here: the most straightforward one implements the conservation of the 
clump's mass 

\ l + 5 n e (r ) ] 1/3 
a = a — — — -— , (19) 
|_ 1 + d n e [r) _ 

where ro is the initial radial distance of the clump from the centre, ao and So are the initial 
clump's radius and overdensity, respectively. On the other hand, as a clump becomes dense 
and moves through the ICM, it may lose its spherical shape, or even break up into several 
smaller clumps. The fragmentation (or mass loss) is clearly inconsistent with the conservation 
law (19), so we shall have to implement the description of this phenomenon at some level. 
The process of fragmentation is quite complex and beyond the scope of the present paper, so 
we take a phenomenological stand, and assume that at each radial distance from the centre 
r the mass of the clump m is reduced with respect to its initial value mo according to the 
prescription 



m 



m (r/r ) fc , (20) 
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where the mass loss index k ^ is a constant. This prescription is equivalent to 

a - a (—} 1 + ne ( r °) 
\r J 1 + 5 n e (r) 

which reduces to the mass conservation (19) if k — 0. Note that the decrease of a implied by 
Equation (21) has the same effect to increase the cross section to volume ratio, which mimics 
the clump's loss of spherical symmetry, with the major axis aligned with the direction of 
motion. 

Note that Equation (20) properly describes the effects of loss of shape and/or break 
up if the clump does not oscillate across the ICM, since in this case m(t) would not be a 
monotonic function of time. 

In Section 5 we shall solve numerically the system of equations (17) separately for the 
prescriptions (19) and (21). In Appendix A we present an analysis of the equilibrium solutions 
of Equations (17) (with the conservation condition 19), and also provide an approximate 
stability criterion for the evolution of a clump (see also Balbus & Soker 1989; Tribble 1991). 



(21) 



4 CHOICE OF THE CLUSTER PROFILES 

As in Paper I, we focus on a specific cluster, i.e. Virgo. We adopt the entropy, temperature 
and electron density profiles derived by Cavagnolo et al. (2009) 2 . We fit the 3D electron 
density with a double beta model 

Ul 11,2 
Ue= [1 + (r/n) 2 ]^ + [1 + (r/r 2 )Y 2 ' 

The 3D entropy index profiles are modelled by Cavagnolo et al. (2009) with the simple 
formula 

K = K + K 1 ( r - ^ . (23) 

These entropy profiles have been deduced from non-deprojected (2D) temperatures: ac- 
cording to Cavagnolo et al. (2009) this does not affect the determination of the 3D entropy 
profiles appreciably. The temperature profiles provided by Cavagnolo et al. (2009) are 2D. 
Since Virgo is almost isothermal (see e.g. Ghizzardi et al. 2004), we do not expect that 
the use of 2D instead of 3D profiles will severely affects our results. We then use the 2D 



2 Their data (for a sample of 239 clusters observed with Chandra) may be found at 
http: //www.pa.msu. edu/astro/MC2/accept/ . 
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2D temperature 


T 
( keV) 


Ti 
( keV) 


n 

( kpc) 


7 






1.893 


3.022 


29.36 


1.543 






3D electron density 


"1 

(10 -2 cm" 3 ) 


ri 
(kpc) 


01 


™2 

(10~ 2 cm" 3 ) 


( kpc) 


02 


1.320 


8.762 


0.24 


15.08 


1.896 


0.819 


3D entropy index 


K 
( keV cm 2 ) 


i^i 
( keV cm 2 ) 


a 








3.53 


146.6 


0.8 









Table 1. Best fitting parameters for the (projected) temperature, and 3D electronic density and entropy index for Virgo. 
More details in the text. 

temperature profile as an estimate of the real 3D profile, and we fit it with the expression 

T = T ° + ™f, (24) 
l + (r/r,)i 1 ' 

already used by Vikhlinin et al. (2006) to model the temperature of galaxy clusters in their 
cool cores. The profiles (22) and (24) are used to determine the hydrostatic gravitational 
acceleration 

GM(<r) k B T fd\nn e dlxiT 



9 



ji m u r \ dlnr dlnr 



(25) 



where M(< r) is the gravitating mass within the radius r, m u is the atomic mass unit 
and fi ~ 0.6 is the mean molecular weight. Table 1 presents the best-fitting values of the 
fit parameters for the electron density, temperature and entropy profiles. The profiles for 
Virgo of T, n e , K and M(< r) are plotted in Fig. 2, compared to the original data from 
Cavagnolo et al. (2009). 

Throughout this paper we assume that the metal abundance is 0.3 times the solar value 
(Cavagnolo et al. 2009), with the heavy element ratios taken from Anders & Grevesse (1989). 



5 NUMERICAL CALCULATIONS 

In this section we solve numerically equations (17) with the temperature, entropy and density 
profiles presented in Section 4. We adopt the embedded Runge-Kutta Prince-Dormand (8,9) 
method (e.g., Press et al. 2007) implemented in the GNU scientific library 3 . The cooling 
function A has been calculated with the XSPEC package for the APEC atomic tables. 

In the next two subsections we present the evolution of the clumps with and without 
angular momentum. 

3 http : //www. gnu. org/sof tware/gsl/manual/html_node/Ordinary-Dif f erential-Equations .html . 
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r (kpc) r (kpc) 

Figure 2. Best fitting profiles for the electron density (top left), temperature (top right), entropy index (bottom left) and 
gravitating mass (bottom right) profiles adopted for Virgo, superposed to the data by Cavagnolo et al. (2009). 



5.1 Zero Angular Momentum 

We follow the evolution of a clump of initial radius do = 100 pc released from rest at 
a distance r = 20 kpc from the cluster's centre. We distinguish the case where the clump 
conserves its mass during the inflow, from the case where the clump undergoes fragmentation. 

5.1.1 Constant Mass Clumps 

The size of a constant mass clump is given, instant by instant, by the conservation con- 
dition (19). We may distinguish three regimes for the accretion of such a clump: stable, 
relatively unstable and absolutely unstable. 

These regimes are exemplified in Fig. 3, showing the evolution of a clump for different 
initial values of the overdensity: 6q = 1.4 (solid line), 5q = 2.0 (dashed line) and 5q = 3.0 
(dotted line). The top left panel of Fig. 3 shows the evolution of the overdensity defined 
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Figure 3. Tracks of a clump with an initial radius a = 100 pc released from ro = 20 kpc from the centre. The lines refer 
to different initial overdensitics: <5o = 1-4 (solid line), <5o = 2.0 (dashed line) and So = 3.0 (dotted line). Top left: overdensity 
vs. time; top right: radial distance vs. time; bottom left: fall velocity vs. time; bottom right: Mach number of the radial fall 
velocity vs. overdensity. 



by Equation (7). The clump with 5 = 1.4 is stable. Fig. 3 shows that its overdensity 
soon drops to zero, oscillating around zero with the Brunt-Vaisala frequency. The clump 
eventually reaches an equilibrium position at r ~ 0.3 kpc. The trajectory of this clump in 
the (5 — A4) plane (where M. is the ratio of the radial velocity to the local sound speed, i.e. 
the Mach number) winds up around the origin, which is eventually reached when the energy 
of the oscillations is dissipated by the drag force. The case of 5$ = 2.0 is an example of a 
relative instability: the overdensity decreases, but the clump is eventually accreted, since it 
reaches the centre of the cluster before being thermally stabilised. If the overdensity exceeds 
a critical value 5c — 2.4 (slightly less than the estimate 5c = 1.7 calculated in Appendix A) 
the overdensity always increases, and the clump is eventually accreted to the centre: this is 
the absolute instability, shown in Fig. 3 by the case 5$ = 3.0. 
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The isobaric prescription for the evolution of very dense (i.e., absolutely unstable) clumps 
breaks down when the sound crossing time of the clumps becomes longer than its cooling time 
r sound > r cooi (Burkert & Lin 2000). For our clump with <5 = 3.0 this happens when 5 ~ 32, 
or T' ~ 0.06 keV. The growth of the overdensity stalls after the equivalence t' CO(A ~ Tg 0und 
has been reached, and from this time on the clump evolves isochorically. The subsequent 
evolution of the clump's temperature is difficult to estimate, but we may note that the cooling 
function sharply drops below T ~ 10 4 K (see e.g. Spitzer 1978; Sutherland & Dopita 1993), 
and since r' cooX greatly increases, the clump is unlikely to cool much below this limit before 
being eventually accreted. The conclusion is that the temperature of the clumps feeding the 
AGN may be T' ~ 10 4 K, but probably not much colder. 

On the other hand, a very dense clump may undergo a gravitational instability: it breaks 
up, decouples from the ICM and eventually freely falls to the centre. This evolution is rapid, 
and it may be possible that the gas of the clump encounters an atomic or molecular phase, 
as suggested by the presence of a large amount of molecular clouds far from the centre of 
clusters (Salome et al. 2008). In this case, it is likely that the gas will cool much below 
1000 K, possibly down to few tens K. 

We have neglected the feedback activity of the AGN, that heats the ICM, but not the 
clumps. This activity increases the ICM entropy, as well as the entropy contrast between 
the clump and the ICM. This effect makes the clumps more unstable than found here. 

The results obtained thus far are in good agreement with those presented in Paper I. The 
small differences are due to our different choices of the metal abundance and the cluster's 
temperature, density and entropy profiles. 

5.1.2 Effects of Fragmentation 

We include the description of the clumps' break up by setting the mass loss index k (Equa- 
tions 20 and 21) to a positive value. Fig. 4 compares the evolution of 5 for clumps with 
different k, but the same initial radius (do = 100 pc) and overdensity (Sq = 1.4); all the 
clumps are released from rest at the distance ro = 20 kpc from the centre. The curves in 
the plot refer to the mass loss indexes for the values k — 0, k — 1.0, k = 2.0 and k = 4.0. 
The most relevant effect of a non- vanishing A; is a strong increase of the relative drag force, 
which suppresses the oscillations for k > 2. 

Fig. 5 presents the evolution of a clump with Sq = 1.4 and mass loss index k = 2.0. For 
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Figure 4. Evolution of the overdensity (loft panel) and the mass m of a clump with respect to its initial value mo (right 
panel) for different values of the mass loss index k (eq. 20): a constant mass clump (k = 0.0; solid line), k = 1.0 (dotted line) 
k = 2.0 (dashed line) and k = 4.0 (large fragmentation of clump; dash-dotted line). The clump has the initial density contrast 
So = 1.4, radius ao = 100 pc and is released with zero velocity and angular momentum from ro = 20 kpc. 

comparison, we have also drawn the tracks of a clump with the same parameters, but k = 0. 
If k = 2.0 the drag force suppresses the oscillations and limits the fall velocity. The clump 
falls monotonically to the centre, and the initial decrease of its overdensity is reversed. We 
do not present the plots for denser clumps since the general pattern remains the same: the 
stronger drag force accelerates the evolution of the clump somewhat, with no other relevant 
qualitative effects. 

Fig. 6 shows the effect of the size a of the clump on its evolution. Since the drag force 
behaves as Fjj oc a~ l , it effectively damps the oscillations of small clumps, whose distance 
from the centre then decreases monotonically. The ratio between the drag force and the force 
of gravity (for the same density contrast) is larger for smaller clumps. Therefore, the trend 
of decreasing clump's radius is like that of breaking (k > 0). Smaller clumps fall slower, and 
therefore their density contrast has time to grow, making these clumps less stable than the 
larger ones. On the other hand, we also expect that thermal conduction may play a role in the 
evolution of small clumps. As discussed in Paper I, we expect that the clumps are protected 
against the evaporation by magnetic fields, which suppress the thermal conduction with the 
ICM by a factor ~ 10~ 3 (Nipoti & Binney 2004) with respect to its nominal value (as given 
by Braginskii 1965 or Spitzer 1962). Magnetic shielding, however, becomes inefficient for 
very small clumps, that will then be evaporated. 

According to our phenomenological recipe (20), for moderate values of k the clumps may 
reach the distance r ~ 1 kpc still retaining an appreciable fraction of their initial mass. 
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Figure 5. Dashed line: tracks of a clump with mass loss index k = 2.0. The initial radius and overdensity ao = 100 pc, <5o = 1.4. 
The clump is released from ro = 20 kpc with zero initial velocity and angular momentum. Top left: overdensity vs. time; top 
right: radial distance vs. time; bottom left: fall velocity vs. time; bottom right: fall Mach number vs. overdensity. The solid 
lines show the evolution of a clump with the same initial parameters, but with k = 0. The case k > favours the accretion 
process. 



Although a more refined analysis is required to make reliable quantitative predictions, we 
are confident that this result is valid in general. At smaller radii (below r ~ 1 kpc) the 
evolution of the clumps is different from the scenario sketched so far, as discussed later in 
Section 6. 



5.2 Clumps with Angular Momentum 

We now consider the accretion of clumps endowed with an initial orbital angular momentum. 
We solve the full set of equations (17), and set the initial angular momentum to 

So 



,1/2 



9o rl 



l + So 



(26) 
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Figure 6. Evolution of ovcrdensity (left panel) and distance from the centre (right panel) of clumps of different sizes and mass 
loss indices. The solid and the dashed line show the evolution of clumps with ao = 20 pc and ag = 100 pc, respectively. For 
both these clumps it is k = 0. For comparison, the dotted line shows the evolution of a clump with k = 2.0 and ao = 100 pc. 
In all cases the initial overdensity is 5q = 1.4, and the clump is released from rg = 20 kpc with zero initial velocity and angular 
momentum. Smaller clumps fall slower, and therefore their density contrast has time to increase. They are less stable and more 
prone to the accretion process. 



for which a clump with zero initial radial velocity also has zero radial acceleration (see 
equation 17c, where go = g(r )). The clumps have the same parameters as in Section 5.1, 
i.e. ao = 100 pc, r = 20 kpc and overdensities Sq = 1.4, 2.0 and 3.0. A summary of the 
evolution of some relevant parameters is plotted in Fig. 7. The evolution of the overdensity is 
not significantly different from that found with zero angular momentum for all the clumps. 
As seen from the shape of the orbit (bottom left panel of Fig. 7) and from the angular 
momentum loss time (Fig. 8) all the clumps quickly lose their angular momentum, and their 
orbit becomes a straight line pointing to the centre. 

The clump with 5 = 1.4 is still stable: it reaches an equilibrium radius r cq ~ 0.3 kpc and 
oscillates around it. The denser clumps quickly become very overdense and are eventually 
accreted by the central black hole. The last panel of Figure 7 shows that the clumps loose 
their angular momentum before cooling. For even denser clumps this may not be the case. 
Fig. 9 presents the case of a clump with initial overdensity 5q = 5. After an initial phase 
during which also this clumps loses angular momentum, the orbit starts to wrap around the 
centre (x, y) = (0,0), with decreasing average distance. This clump, unlike the less dense 
clumps, cools down before losing its angular momentum, so its overdensity 5 soon grows 
very large. The drag force (proportional to (1 + 5)~ 2 ^ 3 ) becomes very small, increasing 77 
and delaying the eventual accretion to the centre. Albeit formally correct, this solution may 
be unphysical for the following reasons. First, the orbital periastron of such dense clumps 
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Figure 7. Tracks of a clump with initial radius ao = 100 pc released from ro = 20 kpc from the centre and having an initial 
angular momentum according to Equation (26). The lines refer to different initial ovcrdcnsities: So = 1.4 (solid line). 8q = 2.0 
(dashed line) and So = 3.0 (dotted line). Top left: overdensity vs. time; top right: radial distance vs. time; bottom left: shape 
of the clumps' trajectories in the orbital plane (the cluster centre is at (x,y) = (0,0)); bottom right: evolution of the ratio 
between the clump's cooling time and the angular momentum loss time. 



is small (below 100 pc). As explained in Section 6, at these radii the collisions with other 
clumps cause the clump to efficiently lose its orbital angular momentum. This accelerates the 
accretion process as required by the cold feedback mechanism. Second, when the overdensity 
value of ~ 45 is reached (for this example, this occurs at t ~ 4.1 x 10 7 yr, when the clump 
is still r ~ 16 kpc away from the centre), the cooling becomes isochoric rather than isobaric 
(Burkert & Lin 2000). The clump's compression lags behind cooling and leaves the size of the 
clumps larger, increasing somewhat the drag force relative to what our calculations assume 
at this phase. 

We may summarise the main results of this section as follows: the angular momentum 
is quickly lost by many clumps, and their behaviour is not much different from the case 
with zero angular momentum presented in Paper I and in Section 5.1. This result removes 
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Figure 8. Evolution of the angular momentum of a clump with initial radius ao = 100 pc released from tq = 20 kpc from the 
centre and having initial angular momentum as in Equation (26). Left panel: the evolution of the angular momentum loss time 
(given by Equation 15); right panel: ratio of the angular momentum l(t) with respect to its initial value 1(0). The lines refer 
to different initial overdensities: So = 1.4 (solid line), <5o = 2.0 (dotted line). So = 3.0 (dashed line) and <5o = 4.0 (dot-dashed 
line). The spikes in the curve Sq = 4 are due to the modulation of the clump's velocity on an eccentric orbit. 




x (kpc) t (yr) 

Figure 9. Evolution of a dense (<5q = 5) clump with initial radius ao = 100 pc released from ro = 20 kpc from the centre and 
having initial angular momentum as in Equation (26). Left panel: the trajectory in the orbital plane. Right panel: ratio of the 
cooling and angular momentum loss time (solid line); ratio of the cooling and sound crossing time (dashed line). 

the main obstacle that stood in the way of the cold feedback mechanism. Namely, we have 
shown that the response of the AGN accretion to cooling at large distances is as rapid as 
required. 



6 COLLISIONAL BRAKING 

In this section we estimate the role of collisions in the innermost region of the cluster. As 
the clumps fall in, they are likely to break up into smaller clumps and/or lose their spherical 
shape. Clumps that break up into very small pieces are likely to be evaporated by thermal 
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conduction (Nipoti & Binney 2004). On the other hand, clumps larger than ~ 10 — 100 pc 
(for a thermal conduction suppression factor 10~ 2 — 10 -3 , Fig. 3 of Nipoti & Binney 2004) 
will survive, and their subsequent evolution is likely dominated by mutual collisions, as 
explained in this section. Loss of spherical shape and fragmentation increase the cross section 
to volume ratio of each clump. With the higher clumps' number density near the centre, 
collisions between clumps are bound to occur. 

Consider the clumps within a sphere S of radius r centred on the AGN. The collision 
time between clumps here is 

r coll = ° ^rei n c (27) 

where v m \ is the relative velocity of two clumps, n c the number of clumps per unit volume 
within S, and a is the collisional cross section. We take it to be the geometric cross section: 
a = 7r(2a) 2 , where the extra factor 2 comes from the finite size of both the target and the 
bullet clump. This time is to be compared with the clump's crossing time of S 

Across ~ r/v c , (28) 

where v c is a typical clump's velocity. Random collisions are frequent if r co n < r cross , and may 
efficiently cause the clumps within S to lose their angular momentum. The ratio r co u/r cross 
therefore estimates the efficiency of the collisional loss of angular momentum; it may be 
expressed in terms of the clumps' volume filling factor ey = 4ira 3 n c /3 as 

Tcoii v c a 1 a/r 

~ ~ — - — , (29) 

r cross v re \ r 3 €y 3 €y 

where we have assumed v c ~ v Te \ in the last equality. For example, if ~ 10% of the clumps 
lose their angular momentum by collisions, then this ratio should be ~ 10. For r/a = 300 
(say, r = 30 kpc and a = 100 pc, or r = 3 kpc and a = 10 pc) the required efficiency 
demands ey ~ 10~ 4 . Namely, it is sufficient that only 0.01% of the inner volume is occupied 
by dense clumps in the inner region. In practice, we expect ey to be larger, e.g. 10~ 3 — 10~ 2 , 
depending on the mass accretion rate. 

We conclude that even a small filling factor is enough to make the collisional loss of 
angular momentum an efficient process. The overall picture is that the drag force dominates 
the angular momentum loss at large radii (r > 1 kpc). At smaller radii (r < 1 kpc) collisions 
dominate, causing the clumps to lose their residual orbital angular momentum. 
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7 CONCLUSIONS AND SUMMARY 



We have studied the role of the cold clumps' angular momentum in the framework of the cold 
feedback model, generalizing the results of Paper I that did not include angular momentum. 
We first studied falling clumps with zero angular momentum. The efficiency of accretion 
increases (i.e., clumps are less stable) as the size of the clump decreases (for a given density 
contrast) and its cross section increases (for a given volume and density contrast), as evident 
from Figs. 4-6. The reason is that a relatively large cross section implies slower infall, such 
that the clump has time to increase its overdensity due to radiative cooling. 

We then include angular momentum. As it is apparent from Figs. 7 and 8, at large radii 
the drag force is the main agent cause for clumps to lose their angular momentum, and 
it is efficient enough to remove the centrifugal barrier of accreted clumps. The clumps can 
easily reach distances of r < 100 pc. At this stage, the clumps may be very dense and cool 
(T < 10 4 K), and they are no more in pressure equilibrium with the surrounding ICM. As 
discussed in Section 6, the collisions between clumps become important, and even for small 
volume filling factors (ey ~ 10 -4 — 10 -2 ) the collision time scale is comparable or shorter 
than the infall time. The residual orbital angular momentum is lost, and the clumps are free 
to accrete on the central galaxy and ultimately feed the SMBH sitting at at the clusters' 
centre. This accretion powers the AGN and activate an efficient feedback with the ICM 
cooling. 

In our calculations the ICM was static. However, along the jets axis there is an outflow, 
and clumps will not fall there freely. We expect most of the clumps to be accreted from 
the equatorial plane. The interaction of the jets and rising bubbles with dense clumps will 
be the subject of a future study. We also note that although the AGN is the main heating 
source, the falling clumps do release gravitational energy (Fabian 2003). 

The solution of the angular momentum problem has removed the main theoretical ob- 
jection to the cold feedback mechanism as a viable model for ICM heating in CF clusters. 
Observational support to the cold feedback mechanism (e.g., Wilman et al. 2009) stand now 
on a more solid ground. 
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APPENDIX A: EQUILIBRIUM SOLUTIONS 

In this appendix we the study the equilibrium solutions of equations (17) for a clump with 
constant mass. At equilibrium the right-hand sides of (17) vanish, which happens if v — 
0, / = 0, 5 = 0; r and <fi are both constant. It is clear that since 5 = the clump is 
indistinguishable from the ICM, so this equilibrium solution describes a clump reabsorbed 
by the ICM. It is worth noting that such an equilibrium solution is allowed only if the cooling 
of the ICM is accounted for, which was not in Paper I. 
The linearised system (17) reads 

%- = u (Ala) 
dt v ; 

du ^ / a ■% i \ 

- = -gS (Alb) 

d5 n d In K u' r , . . 

T77 = f n —7— + 5 ' Alc 

dt 5 dr r coo i 

where 

d In A 



. du 



s=o dlnT 



(A2) 



and 

-1 _ 2 »h n e A{T) 
7-cooi - 5 kBT ( AS ) 

is the inverse of the ICM cooling time. 

The linearised system admits solutions 5 oc u oc e A * if A is the root of the characteristic 

equation 

A f \ 2 - —A + u>ly) = 0, (A4a) 

V T-cool / 

where the (squared) Brunt- Vaisala frequency 

q „ , ^ r -x Q d\n K , L , . 

ul Y = -lg-V\nK=l g -— (A4b) 

is positive in a convectively stable cluster. For typical parameters cubv 3> r C ooi' so the eigen- 



values are 

A = A~cj'/2r coo i± iwbv- (A5) 

One of the eigenvalues is zero, so the equilibrium is non-hyperbolic and a simple lin- 
ear analysis does not suffice to characterise the global stability of the system (see e.g. 
Glendinning 1994). The perturbations oscillate around the equilibrium position with the 
Brunt- Vaisala frequency. If uj' (defined by Equation A2) is positive, the perturbations' amp- 
litude grows on a cooling time. The reason of this (overstable) growth is easily understood: 
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after the clump has completed an oscillation, its density has grown a little on account of 
the radiative losses, which causes the amplitude of its next oscillation to be larger. The 
overstability criterion u' > was already found and discussed by Balbus & Soker (1989): 
in particular see their equation (2.1), with Jzf = n 2 A(T). The oscillations' amplitude of an 
overstable clump will soon grow so large that it is not possible to disregard the non-linear 
terms of the equation, invalidating this simple linear analysis. 

The global behaviour of the system may be studied in the "reduced" set of equations 
made by the velocity and overdensity equations only. In this system r is merely a constant 
parameter. With a suitable scaling of the variables, the reduced system reads 

-j- = u + ^ u(6, T) (A6a) 



1 + 6 dt 
d 

~dt 1 \ l + 6 



1 /3 

l + 5)^ = - V (±±f) ' \u\u-5, (A6b) 



where 



V= 'i C D — £ = (W B V Tcool) 1 . (A6c) 

This system admits the trivial solution 5 = u = 0: this is the only proper equilibrium 
solution of the total system, since the velocity must vanish there. The non-trivial equilibrium 
solutions are 

r / 1 _|_ r \ 1/3 

\u\u= — (-—4) u = -Zu(5,T). (A7) 



1/2 

(A8) 



7] \ l + 6 

The first equation gives the negative (inwards) velocity 



Ui 



6/1 + 6, 



1/3 



T) V 1 + 6 

This is the "terminal" velocity of a clump, resulting from the equilibrium between the 
gravitational pull and the drag force. The overdensity corresponding to the velocity u t is 
implicitly defined by the equation 



\l + <5</ 

where x — V^, 2 an d u t = &(6t, T). This formula will be useful in a moment to derive a simple 
stability criterion for the thermal stability of a clump. 

The study of the eigenvalues of the jacobian matrix to assess the stability of the equilib- 
rium points O = (0, 0) and A = (6f, u t ) is not particularly illuminating. It is more convenient 
to plot (Fig. Al) the directions of the vector fields d5/dt and du/dt given by (A6). The lines 
defined by the equalities d6/dt = (dashed line) and du/dt = (solid line) divide the 
plane (6 — u) into five regions, labelled I, II, III, IV, V. A point starting from the upper 
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Figure Al. Vector fields defined by the reduced equation system (A6). The dimensionlcss parameters (defined in Appendix A) 
are r\ = 1.27 X 10 2 , £ = 1.64 X 10 — 2 and x = 3.42 X 10 -2 . They correspond to a clump of radius ao = 100 pc released from 
rest at the radius tq = 20 kpc in a cluster like Virgo. The points for which dS/dt = lie on the solid line, the points for which 
du/dt = lie on the dashed line. 

right of region I drifts to the lower right of the plane, i.e. its overdensity increases without 
limits, and its velocity approaches the terminal velocity (A8). If the initial position is on the 
left of region I, the trajectory crosses vertically the solid line, moves leftwards to region II, 
and crosses horizontally the dashed line, moving to region IV. The trajectory coasts the 
dashed line, moving to region V and finally to region I again. The overdensity is now smaller 
than the initial one, and this trajectory winds up around the point O, which is stable. This 
behaviour shows the existence of a critical initial overdensity 8c- If 8q > 8c the clump is 
(absolutely) unstable: its overdensity grows more and more, and it is eventually accreted on 
the cluster centre. 

We now derive simple formula for the critical density 8c of a clump released from rest. 
Consider the trajectory of a clump in the (5 — u) plane. The starting point of a clump 
released with zero velocity is S = (5q,0). In a neighbourhood of S the trajectory may be 
approximated by the straight line 

u = u'(8-8 ), (A10) 
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where 

u> = (All) 

is the derivative du/d5 evaluated at S, and ujq = u(5o,T). The line joining S and A has 
(negative) angular coefficient 

m = (A12) 

Of. - o 

where u t = — £ oj t , according to the second equation A7. If w' = m the trajectory originating 
in S hits the point A: the value of 5q for which this occurs corresponds to the critical 
overdensity Sc', if \u'\ < \m\ the trajectory entirely lies in the region I of the (5 — u) plane 
(Fig. Al), and if \u'\ < \m\ the clump moves to region II, and its overdensity decreases. The 
equation v! = m reads 

(l + dcY£uc S t - de- 
Equations (A9) (with S = 5c) and (A13) provide a closed system for 5c and 5 t . This system 

is quite complex, but it may be simplified under the assumption ( < 1, which is true for 

typical parameters. We may therefore expand 

Sc = S t + 1 6$ + e 4 2) + o(0 2 (A14) 
To zero order in £ Equation (A13) simply gives Sc = St, and Equation (A9) becomes St = 
X^t- This equation in St may be solved numerically, with the temperature T entering as a 
parameter. Plugging the expansion (A14) into Equation (A13) and comparing the coefficients 
of the powers of £ we find 5^ = and 8q = bj\ (1 + 5 t ) 2 /5 t , to give 

5c = 5 t + e U ^\ +5t)2 +o(0 2 - (A15) 
5 t 

Since £ ~ 10~ 2 1, the zero-order approximation is usually satisfactory. In this case the 
value 5c — 5 t only depends on the dimensionless parameter x', a couple of examples are 
shown in Fig. A2 for different temperatures. 

Fig. A3 shows the radial profile of Sc — as calculated from the second-order formula (A15)- 
for some values of the clumps' initial radius do- The radial temperature, density and entropy 
profiles are those of the Virgo cluster. 

A final warning about this stability analysis is in order. In the present analysis we have 
neglected the motion of the clump through the cluster: the radius entered only implicitly 
into the coefficients rj and £, treated as constants. If the radial motion of the clump is 
considered, the situation is slightly more involved. First, the clump slightly drifts inwards, 
and its evolution occurs at a smaller distance to the centre, where 5c is larger (if r > 0.2 kpc, 
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Figure A2. Critical overdensity <5t, above which constant mass clumps are thermally unstable. St is plotted against the 
dimcnsionlcss stability parameter \ denned by Equation (A9). The function u) (defined by Equation 13) has been evaluated at 
T = 2.0 keV (solid line) and T = 3.0 keV (dashed line). 
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Figure A3. Critical overdensity <5c, above which constant mass clumps are thermally unstable as a function of the distance 
from the centre, for Virgo. The curves refer to different initial radii of the clumps: ao = 25 pc (solid line); ao = 50 pc (dashed 
line) and ao = 100 pc (dotted line). Clumps with S > 8q are accreted, the other are thermally stabilised at r > 0. 
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Fig. A3). The value of 5c found with the analysis of the reduced system (A6) is therefore 
slightly underestimated. Second, even a clump with 5 < 5c may be accreted to the centre: 
we term this behaviour as "relative instability". If 5q is not much smaller than 5c, the 
clump's overdensity decreases, but the clump may nevertheless reach the centre before the 
equilibrium value 5 = is attained. The value 5c calculated here may only distinguish 
between the "relative instability" (for which 5 < 0, with possible eventual accretion) and the 
"absolute instability" , for which 5 > all the way. 
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